Inducing an Optical Feshbach Resonance via Stimulated Raman Coupling 
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We demonstrate a novel method of inducing an optical Feshbach resonance based on a coherent free-bound 
stimulated Raman transition. In our experiment atoms in a 87 Rb Bose-Einstein condensate are exposed to 
two phase-locked Raman laser beams which couple pairs of colliding atoms to a molecular ground state. By 
controlling the power and relative detuning of the two laser beams, we can change the atomic scattering length 
considerably. The dependence of scattering length on these parameters is studied experimentally and modelled 
theoretically. 



PACS numbers: 34.50Rk, 32.80.Pj, 03.75.Nt, 34.20.Cf 
I. INTRODUCTION 

Feshbach resonances have become a central tool in the 
physics of ultra-cold quantum gases during the last years be- 
cause they allow for a tuning of the interactions between 
atoms. Controlling interparticle interactions is a central key 
in many fields of modern physics and is especially relevant 
for future applications in quantum computation and exploring 
novel many-particle quantum effects. Beautiful experiments 
using magnetically tunable Feshbach resonances 111 lai have 
been performed, ranging from ultra-high resolution molecu- 
lar spectroscopy |3] to the coherent coupling of atomic and 
molecular states jj] as well as the creation of bright matter 
wave solitons |3l ■ It also led to the production of new atomic 
0] and molecular Q] Bose-Einstein condensates (BEC) and 
allowed to control pairing in ultra-cold fermionic gases |8|. 

Recently we demonstrated how atom-atom interactions in a 
87 Rb BEC can also be tuned with an optically induced Fesh- 
bach resonance 0] (see also fioln . a scheme which was orig- 
inally proposed by Fedichev et al. fill fr3l . Optically in- 
duced Feshbach resonances offer advantages over magneti- 
cally tuned Feshbach resonances since intensity and detuning 
of optical fields can be rapidly changed. Furthermore complex 
spatial intensity distributions can be easily produced and op- 
tical transitions are always available even when no magnetic 
Feshbach resonances exist. A disadvantage of optically in- 
duced Feshbach resonance is the inherent loss of atoms due to 
excitation and spontaneous decay of the molecular state |3l ■ 
Typical lifetimes for excited molecular states are on the order 
of 10ns which corresponds to a linewidth of 2n x 16MHz. 
Evidently, coupling to molecular states with longer lifetime 
should improve the situation. Ground state molecules are sta- 
ble against radiative decay, and narrow transition line widths 
on the order of kHz have been observed in two-photon Raman 
photoassociation fl3lfl4ll . This raises the question whether it 
is possible to create optical Feshbach resonances using stim- 
ulated Raman transitions and whether this scheme might be 
advantageous compared to the one photon optical Feshbach 
resonance. 

In this paper we indeed demonstrate that optical Feshbach 
resonances can be induced using a coherent two-color Raman 
transition to a highly vibrationally excited molecular ground 



state in a Rb BEC. In the experiment we show how the scat- 
tering length and loss rates can be tuned as a function of the 
intensity of the lasers and their detuning from molecular lines. 
We use Bragg spectroscopy Jl5ll as a fast method to measure 
the scattering length in our sample To fit and analyze 
our data we use a model by Bohn and Julienne fl^l . We find 
that using the Raman scheme for optically induced Feshbach 
resonances leads to similar results in tuning of the scattering 
length as for the single photon Feshbach scheme. The Raman 
scheme does not lead to an improvement compared to the one- 
photon scheme because its atomic loss rate is not lower for 
a given change in scattering length. However, using a stim- 
ulated Raman transition does offer experimental advantages. 
To tune over the Feshbach resonance, the relative frequency of 
the two laser beams only has to be changed typically by sev- 
eral MHz which can be conveniently done using an acousto- 
optic modulator. This allows for very fast and precise control 
of the scattering length. On the other hand, working with a 
one-photon optical Feshbach resonance in the low loss regime 
typically requires large detunings and scan ranges on the or- 
der of GHz. The Raman scheme relaxes the necessity for ab- 
solute frequency control of the lasers which can be tedious to 
maintain far away from atomic lines. Since off-resonant light 
fields in general lead to dipole forces acting on the atoms, a 
variation of the scattering length via optical tuning leads to a 
variation of the dipole forces on the atomic sample. This un- 
wanted effect can be made negligible for the Raman scheme 
which tunes over resonance within a small frequency range. 

The paper is organized as follows: We start in section ITT1 
by discussing the Raman scheme with a simple theoretical 
model. In section [Hi] we describe in detail our experimental 
setup and the measurement method. In section HV1 we discuss 
the experimental results which are compared with a theoret- 
ical model. The appendix gives details of the model that is 
used to describe the data. 



II. RAMAN SCHEME FOR OPTICAL FESHBACH 
TUNING 

Before discussing optical Feshbach tuning based on a two- 
photon Raman transition, it is instructive to briefly recall the 
one -photon scheme first l9l Ull [l2ll . This configuration uses 
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FIG. 1 : Schematic diagram of the transitions used for optically cou- 
pling the collisional state |0) to molecular states |1) and |2). |1) is 
electronically excited whereas |2) is in the electronic ground state. 
Ai and A2 are defined to be positive for the shown configuration. 
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FIG. 2: A"; ne i and scattering length a according to Eqs. Q and 
plotted for two values of Ai. Solid line: A\/2k = 100 MHz, dashed 
line: h\/2% = 50MHz. The other parameters are Tqi/2k = 50kHz, 
£2l2/2;r = 30MHz, Yi/2z = 25MHz. The wavenumber fej = 2.5 x 
10 corresponds to the finite size of the condensate wavefunc- 
tion. ao is the Bohr radius. 



a single laser beam tuned close to a transition from the scat- 
tering state of colliding atoms to a bound level in an excited 
molecular potential (states |0) and |1) in Fig. 0. Varying 
the detuning Ai or the intensity I\ modifies the coupling and 
hence the scattering length. Atomic loss can occur through 
population of the electronically excited molecular state which 
has a decay width of j\ . 

Introducing a second laser as shown in Fig^will now cou- 
ple the collisional state |0) to a bound level |2) in the ground 
state potential. As we will show, this allows for a tuning of 
the scattering length similar to the one-photon scheme. We 
now have, however, four parameters which can be used to in- 
fluence the scattering length: the intensities I\ and h of lasers 
1 and 2 and the detunings Ai and A2 as shown in Fig. ^1 23]. 

From Eqs. (4.8)-(4.11)] one can extract approximate 
expressions for the inelastic collision rate coefficient K\ nf \ and 
the scattering length a in a Bose-Einstein condensate 12411 : 
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Here Toi denotes the on-resonance stimulated transition rate 
from |0) to |1) and is proportional to l\. £2 12 is the Rabi fre- 
quency for the coupling of the states 1 1) and |2) and is propor- 
tional to -y/TJ. Kk\ is the relative momentum of the collision, 
where Ti is Planck's constant divided by 2k. «b g is the back- 
ground scattering length and m is the atomic mass. 

Equations Q and 10 neglect spontaneous decay from 
state |2) (y 2 = 0) and assume r i <C Yi- Setting i2 [2 = 
yields the expressions for the one-photon Feshbach resonance 
as given in j9}|. Equations Q and (0 yield a Lorentzian and 
a corresponding dispersive line shape as a function of A\. In 



our experiments, however, we hold Ai constant and scan A 2 . 
Figure [2] shows typical curves for Ki ne \ and a for two detun- 
ings Ai. The curves for K me \ are slightly asymmetric, but for 
Ai Q.\2 they can be well approximated by Lorentzians. This 
can be seen by expanding the denominator of Eq. Q in terms 
of A2 at the resonance position. A light shift displaces the po- 
sition of the resonance to Qi 2 /Ai. It is also interesting to 
note that the resonance width decreases with increasing de- 
tuning Aj as yi (^12/Ai) 2 . 

In a sense the two-photon Raman-Feshbach resonance can 
be coined in terms of a one-photon Feshbach scheme. The 
detuning A2 effectively replaces the detuning Ai of the one- 
photon Feshbach scheme 12511 . 

Since Eqs. Q and (|2) have exactly the same form as for 
the one-photon Feshbach resonance, it follows that, given 
a fixed free-bound transition rate Toi, the maximum tuning 
range of the scattering length for the two-photon case cannot 
be larger than in a one-photon scheme. Furthermore, given a 
fixed change-in-scattering-length, the loss rate as determined 
by Ki ne \ is not lower for the Raman scheme than for the one- 
photon scheme. 



III. EXPERIMENTAL SETUP AND METHODS 

A. Production of BEC 

For the experiments we produce 87 Rb BECs of typically 
1.2 x 10 6 atoms in the spin state \F = j,mp = —1). Our setup 
comprises a magnetic transfer line 11711 to transport atoms 
from a magneto-optic trap (MOT) chamber to a glass cell 
where the BEC is produced and all experiments are carried 
out. In a first step about 3 x 10 9 atoms are loaded within 4 s 
into a MOT directly from the background gas and are then 
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cooled further to about 50 pK in a molasses cooling phase. 
After optically pumping into the \F = l,mp = —1) state we 
load the atom cloud into a magnetic quadrupole trap with a 
gradient of 130 G/cm in the (strong) vertical direction. Within 
1 .4 s the atoms are then moved via a magnetic transfer line 
12^1 over a distance of 48 cm including a 120° corner into a 
glass cell which is at a pressure below 10~ u mbar. In this 
cell we finally load the cloud into a QUIC trap HH ending up 
with typically 4 x 10 8 atoms at a temperature of about 250 |jK. 
All three coils of the QUIC trap are operated at a current of 
40 A, dissipating 350 W. This results in trap frequencies of 
ft>radiai/2?r = 150Hz and © ax j a i/27T = 15 Hz at a magnetic bias 
field of 2 G. To achieve Bose-Einstein condensation we use 
forced radio-frequency evaporation for a period of 20 s. The 
stop frequency is chosen so that we end up with condensates 
with a thermal background of about 25% of non-condensed 
atoms. At this value we concurrently get the highest number 
of atoms in the condensate and good reproducibility. For our 
measurements we consider only the condensed atoms. 
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B. Raman lasers 

To realize the Raman scheme shown in Fig.^ we use the 
electronically excited molecular state |1) = |0~,V = l,/ = 2) 
located 26.8 cm -1 below the (S +-P3/2) dissociation asymp- 
tote H dl. About 290 MHz below the 7 = 2 line, there is 
another rotational level with 7 = i27tl . Although about five 
times weaker than the 7 = 2 line, its effect cannot be totally 
neglected in our experiment. We choose level |2) to be the 
second to last bound state in the grou nd state potential. It has 
a binding energy of 636 MHz x h 11311 where h is Planck's con- 
stant. 

The Raman laser beams are derived from a TkSapphire 
laser using an acousto-optical modulator at a center frequency 
of about 318 MHz in a double-pass configuration. This allows 
precise control of their relative frequency difference over sev- 
eral tens of MHz. Both Raman lasers propagate collinearly 
and are aligned along the weak axis of the magnetic trap (see 
Fig. They have a 1/e 2 waist of 76 pm, and their linear 
polarization is perpendicular to the magnetic bias field of the 
trap. 

The TkSapphire laser is intensity stabilized and its fre- 
quency has a line width of about 3 MHz. In order to stabilize 
its frequency relative to the photoassociation lines, the laser is 
offset locked relative to the D2 line of atomic rubidium with 
the help of a scanning optical cavity. This yields an absolute 
frequency accuracy of better than 10 MHz. In all our experi- 
ments the Raman laser intensities were set to I\ — 300 W/cm 2 
and I2 = 60 W/cm 2 at the location of the condensate, if not 
stated otherwise. 



C. Bragg spectroscopy 

To measure optically induced changes in the scattering 
length a, we use Bragg spectroscopy 191 fl5Tl . This method 
allows for a fast measurement on time scales below 100 [is 



FIG. 3: Top left: Experimental arrangement of the laser beams (top 
view). Top right: Absorption image obtained after Bragg-diffracting 
a portion of the atoms to a state with a momentum of two photon 
recoils (lower atom cloud) and subsequent time of flight expansion. 
Bottom: Bragg resonance curves for two different relative detunings 
of the Raman lasers. The relative shift of 700 Hz is due to two differ- 
ent scattering lengths which are optically induced in the condensates. 
The atom numbers are the same for both curves. Shown is the per- 
centage of diffracted atoms versus the frequency difference of the 
Bragg lattice beams. For better comparison we have scaled up the 
right curve by 10 %. 

which is vital because of the rapid photoassociation losses we 
experience in our experiments. A moving lattice composed 
of two counter propagating beams with wavenumber k and a 
frequency difference Af is used to diffract some of the con- 
densate atoms to a state with non-zero momentum. When 
energy and momentum conservation are fulfilled, the Bragg 
lattice resonantly transfers a momentum of two photon recoils 
2hk in a first order diffraction process. For the case of a ho- 
mogenous condensate of density n, the resonance energy for 
Bragg diffraction is given by the sum of transferred kinetic en- 
ergy hAfo = (2hk) 2 /2m and the change in mean-field energy 
Anh 2 na/m 12811 . This corresponds to a frequency difference 
of the Bragg lasers of 

2h 

A/ r = A/ H na. (3) 

m 

If the condensate is initially not at rest, the kinetic energy 
contribution A/o to the Bragg resonance frequency (Eq. [3ji 
contains an additional term 2hkp/m, where p is the initial 
atom momentum in the direction of the Bragg lattice. In our 
experiments we observe such a motional shift corresponding 
to condensate momentum of up to p = 0.1 hk. This momen- 
tum can partly be attributed to optical dipole forces of Ra- 
man beams which are slightly non-centered on the conden- 
sate. Partly it can be attributed to a forced oscillation of the 
condensate in the magnetic trap at 150 Hz which coincides 
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with the trapping frequency. Since this oscillation is driven by 
a higher harmonic of the line frequency (50 Hz), it is in phase 
with the line frequency and we are able to stabilize the ini- 
tial condensate momentum by synchronizing the experiment 
to the line. A stable initial condensate momentum can then 
be determined and canceled out by measuring A/ r alternately 
for Bragg diffraction to the +2hk and —2hk momentum com- 
ponent. After these measures we were left with a residual 
momentum noise level of up to p = 0.01 Tik. 

In our setup the Bragg lattice beams are oriented along the 
horizontal direction perpendicular to the Raman laser beams 
(see Fig. [3} and have a width of « 0.9 mm. We extract both 
beams from a single grating stabilized diode laser and use 
two acousto-optical modulators to control the frequency dif- 
ference. The laser is tuned 1.4nm below the 87 Rb line 
which defines A/o in Eq. 0to be 15. 14kHz. This frequency is 
much larger than the typical mean field contribution, 2hna / m, 
which in our experiments was below 3 kHz. 

We illuminate the trapped condensate for 100 |js with 
the Bragg lattice light. After 12 ms of free expansion the 
diffracted atoms are spatially separated from the remaining 
atoms. Absorption imaging allows us to determine the diffrac- 
tion efficiency. By adjusting the Bragg laser intensity (typi- 
cally 1 mW) we keep the maximum diffraction efficiency be- 
tween 15 % and 20 %. When we scan the frequency difference 
Af and measure the fraction of Bragg diffracted atoms we ob- 
tain curves as shown in Fig. [5] (bottom). These curves have 
a width of approx. 9 kHz as determined by the 100 |js length 
of our Bragg pulses. The shape of the curves is given by the 
Fourier transform of our square light pulses which we use to 
fit the data to obtain the resonance position A/ r |J. The shift 
between the two Bragg spectroscopy curves in Fig.0(bottom) 
is optically induced by shining in the Raman lasers at the same 
time as the Bragg lattice. For both curves the atom numbers 
are the same and Ai = 60 MHz. Only the Raman detuning A2 
differs by 26 MHz. According to Eq. Q this observed shift 
in Bragg resonance frequency is then due to a change in scat- 
tering length, induced by tuning At. This demonstrates that 
we can tune the scattering length a with a Raman Feshbach 
resonance. 



D. Determination of Scattering Length 

We use Eq. Q to determine the scattering length a from the 
measurements of the Bragg resonance frequency A/ r . Equa- 
tion 0, however, is derived for the case of a homogenous con- 
densate. Our trapped condensate, in contrast, which is subject 
to photoassociation losses exhibits a time and position depen- 
dent density n. This can be taken into account by replacing 
the density n in Eq. Q by an appropriate effective value h. 

A simple approach to estimate n is to calculate the spatial 
and time average of the condensate density n over the dura- 
tion of the Raman pulse length T. For this we use the rate 
equation for the local density n = — 2/fj ne i n 2 for two-atom 
losses. The inelastic collision rate coefficient governing 
this process is obtained from measuring the atom number at 
the beginning and the end of the light pulse. This procedure 



already yields good results which differ less than 10% from 
an improved approach which we use for our data analysis and 
which is explained in the following. 

The improved approach consists of a full numerical simu- 
lation which describes Bragg diffraction in a dynamically and 
spatially resolved way. We divide the condensate into density 
classes and and treat their time dependence individually. The 
Bragg diffraction process is identified as a Rabi oscillation be- 
tween a coherent two level system, i.e., the BEC component at 
rest and the Bragg diffracted component. The changing den- 
sity of the condensate due to loss is reflected in a time depen- 
dent resonance frequency (see Eq. Q)- As a result of these 
calculations we obtain for each density class a Bragg reso- 
nance curve similar to the experimental ones shown in Fig. [3] 
Averaging over these resonance curves and determining the 
center position yields the simulated value for the Bragg reso- 
nance A/r. Using Af r = A/o + 2Ma / m we can then determine 
the effective density h. 



IV. RESULTS 
A. Raman Scans 

Figure presents measurements where the detuning Ai of 
laser 1 from the excited molecular state is set to A\/2n = 
60 MHz. The intensities of the Raman lasers 1 and 2 are 
300W/cm 2 and 60W/cm 2 , respectively. Fig. a) shows 
the atom number after illuminating a condensate of initially 
1.4 x 10 6 atoms for 100 |js with the Raman lasers. Scanning 
the Raman detuning A2 we find a strong loss of atoms on res- 
onance. As already expected from Eq. Q the line shape is 
slightly asymmetric. Figure 0b) shows the resonance fre- 
quency A/ r as measured by Bragg spectroscopy. When we 
analyze the data in Fig. 0a) and (b) with the improved proce- 
dure described in section lTXl Dl we obtain values for the scatter- 
ing length which are shown in Fig. 0c). The scattering length 
a shows a dispersive variation between 50 «o an d 140 ao as 
we scan over the resonance. The dispersive scattering length 
curve is offset by about 20 «o from the background scattering 
length ahe= 100ao f° r 87 Rb in the \F = 1,/n/? = — 1) state 
I20l I21L I22I1 . This is due to the one-photon Feshbach tuning of 
laser 1, in agreement with our previous measurements |9|. 

We find that Eqs. Q and (|2ji are not sufficient to describe 
these data properly, mainly because they neglect the decay rate 
72- A more complete model (see Appendix), also taking into 
account both the / = and 7 = 2 rotational levels, was used 
for creating fit curves ll29il . depicted as solid lines in Fig.0 
The fact that the data for atomic loss as well as for the scatter- 
ing length a are both well described by the theoretical curves 
is an intrinsic consistency check for our model and our data 
analysis. 

The shape of the signal A/ r in Fig. 0b) is a combination of 
the effects of the varying scattering length a and the varying 
atom number (see Eq. (|3}). This is illustrated by the dashed 
and continuous lines in Fig. 0b): The dashed line shows the 
expected signal if only the variations in atom number would 
occur and the scattering length stayed constant ll30ll . The con- 
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FIG. 4: Optical Feshbach resonance using a Raman scheme, (a) 
shows the measured atom number after the Raman pulse, (b) the 
measured Bragg resonance frequency and (c) the scattering length, 
as determined from (a) and (b). In (a) the filled circles correspond to 
measurements where Bragg spectroscopy was used to determine the 
scattering length, while the small open circles stem from additional 
loss measurements without Bragg spectroscopy. From our measure- 
ments we estimate the uncertainty of the Bragg resonance frequency 
to be smaller than ±100 Hz, as indicated by the error bar in (b). The 
solid lines in (a), (b) and (c) are from a model calculation (see Ap- 
pendix and text). The dashed line in (b) shows the expected signal if 
there was only loss in atom number but no change in scattering length 
(see also discussion in text). The vertical line indicates the location 
of maximal loss in (a) and helps to compare the relative positions of 
the three curves. 

tinuous line takes both the variations in atom number and in 
scattering length into account. The deviation of the measured 
data points from the dashed line is due to an optical induced 
change of the scattering length. 

B. Dependence on detuning 

We now investigate how detuning Ai affects the scattering 
length a. Figure |5] shows a set of three curves showing the 
scattering length for detunings Ai /2k = 40, 60 and 90 MHz. 

The measurements clearly show that the position and width 
of the resonances depend on Ai. The change in position can 
be mainly explained as light shifts of levels |1) and |2) due 
to laser 2. The decrease of the resonance width with increas- 
ing detuning A] follows directly our discussion in section ITfl 
The solid lines are model calculations as described in detail 
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FIG. 5: Variation of the scattering length with Raman detuning for 
three various detunings Ai from the excited molecular state. The 
solid line is a calculation (see Appendix) which uses a single set of 
parameters for all curves. 
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FIG. 6: Maximum variation in scattering length Aa = a max — a m ; n 
versus one-photon detuning Ai, Solid line: full model calcula- 
tion (see Appendix). Dotted line: 3-level model (see Fig.0, with 
Yi/2k = 2MHz. Dashed line: 3-level model, with 72/2^ = 100 kHz. 

in the Appendix. They are derived from a simultaneous fit to 
the data shown in Fig.[5]and a large number of atom loss mea- 
surements with different detunings (not shown). The set of fit 
parameters is listed in Appendix IA 21 We also use this same 
set of parameters for the theoretical curves in Figs.[6]and[7] 

It is interesting to note from Fig.yjthat the amplitude of the 
dispersive scattering length signal decreases as Aj becomes 
larger. This is not to be expected from the simple model 
Eqs. and (0. To investigate this effect we have performed 
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scans for atom loss and scattering length for several detun- 
ings Ai. Figure [5] shows the maximum variation in scattering 
length, Aa 

— flmax — ^mini obtained for detunings A\ ranging 
from 40 MHz to 200 MHz. Here, a max and a m i n are the max- 
imal and minimal scattering length values for corresponding 
scan curves. Typical scan curves are shown in Fig. [5] Each 
data point in Fig. [6] was derived from a complete scan and 
corresponds to one day of data collection. 

An analysis of our data using our theoretical model indi- 
cates that the decrease of Aa as a function of Ai is a conse- 
quence of two effects. 

(i) To properly model these measurements we have to as- 
sign to the molecular state |2) in the ground state potential 
a non-negligible decay width j2/2n m 2 MHz. For compari- 
son, two calculations of a 3 -level model are plotted in Fig.|6] 
For small 72/2 7T = 100 kHz (dashed line) Aa decreases only 
weakly. For Ji/ln = 2 MHz (dotted line) the theory fits the 
data much better. Such a large decay rate of a ground state 
level is surprising. It seems too large to be explained purely 
by collisions. We find that the decay rate increases with the 
light intensity. At low light powers of a few W/cm 2 we have 
observed very narrow linewidths 72/2 ft on the order of a few 
kHz, similar to the values reported by 

mill- 

The broad- 
ening of the molecular ground level could be due to coupling 
to excited molecular levels. We can exclude, however, from 
our experimental data that these levels are located within our 
experimental scanning range between the states |1) and |3). 
This would lead to additional resonance features in the scat- 
tering length, absorption and light shifts, which are inconsis- 
tent with our data. In contrast, our data indicate a relatively 
constant background loss rate of the ground level over the ex- 
perimental scan range. This allows us to analyze the data suc- 
cessfully with our simple few-level model. Besides coupling 
to excited molecular states, we suspect that coupling to the 
of-wave shape resonance of the scattering channel also gives 
rise to a sizeable contribution to the molecular decay rate. Be- 
cause the c/-wave shape resonance is located very close (a few 
MHz) to threshold, it is resonantly coupled to the molecular 
ground state level via the Raman transition. To include the 
shape resonance is beyond the reach of our simple model and 
has to be investigated later. 

(ii) The second reason for the decrease in Aa is a quan- 
tum interference effect involving both the J = 2 and / = 
rotational levels as predicted by our model. At a detuning of 
A\/2n « 250MHz the interference effect leads to a complete 
disappearence of the optical Feshbach resonance. We observe 
this in a corresponding disappearence of the atom loss fea- 
ture in our measurements (not shown). The interference effect 
alone, i.e., without a 2 MHz linewidth, is not sufficient to ex- 
plain the experimental data in Fig. [5] 



C. Dependence on intensity 

From the simple model Eq. (|2j it is clear that the maximum 
variation in scattering length Aa is proportional to Tqi and 
consequently scales linearly with the intensity I\ of laser 1 . 
We have verified this dependence recently 0] for the case of 
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FIG. 7: Maximum variation in scattering length Aa = a max — a m ; n 
versus h. For this data set /] = 300 W/cm 2 and Ai = 60 MHz. The 
solid line is a full model calculation (see Appendix). The dashed line 
stems from the same model, but with Yi/2k set to 100 kHz and is 
scaled by a factor of 0.84 for better comparison. 



a one-photon optical Feshbach resonance. 

In contrast, the dependence of Aa on intensity I2 of laser 2 is 
not so trivial. According to the simple model Eqs. Q and Q 
which neglect the decay rate 72, the maximum change Aa is 
independent of h- It is also clear, that for h — we have 
Aa = since there is no dependence of scattering length on 
A2 at all. This unphysical discontinuous behavior can be re- 
solved, if we introduce a finite decay rate 72 > 0. We then 
find that for increasing intensity h, Aa rises continuously from 
zero to a value where it saturates. We observe this general 
behavior in our measurements presented in Fig. Q Our full 
model, as described in the Appendix, describes the measured 
data well if we set the decay rate to 72/2^ = 2 MHz (solid 
line). In contrast, the dashed line in Fig.Qshows the calcula- 
tion for the same model where ji is set to 72/271: = 100 kHz. 
Saturation then occurs at a much lower intensity I2 than for 
7 2 /27T = 2MHz. 



V. CONCLUSION 

Our experiments demonstrate the use of an optical Fesh- 
bach resonance for tuning of the scattering length via stimu- 
lated Raman coupling to a bound molecular state. Our results 
show that there is no advantage over a one-photon scheme 
when comparing tuning range and loss rate. However, for cer- 
tain applications a Raman scheme is experimentally more fa- 
vorable since it demands a lower tuning range of the lasers. 
Our presented theoretical model is in good agreement with 
our data and might be helpful when tailoring experimental pa- 
rameters for a specific application. Furthermore it gives in- 
sight into the process of creating stable ultracold molecules 
via two-photon photoassociation. 
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FIG. 8: Extended level scheme (compare to Fig0 for the 4-level 
model (right hand side) and its extension with 6 levels. State |3) 
corresponds to the 7 = level and lies 290 MHz below the J = 2 level 
|1). The 4-level model is based on levels |0), |1), |2) and |3). The 
auxiliary levels 1 1') and |3') are introduced in the extended model to 
describe the coupling between \2) and |1), |3), respectively, due to 
laser 1 (see text). 
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APPENDIX A: THEORETICAL MODEL AND FIT 
PARAMETERS 



We use a theoretical model by Bohn and Julienne [16] to fit 
the data in Figs. |4jQ In the following we give a short sum- 
mary of this model and present the procedure to calculate the 
scattering matrix S, the loss coefficient K me \ and the scattering 
length a. The model has the advantage that it is concise and 
intuitive and it allows treatment of multilevel systems with 
several couplings between the levels. The numerical calcula- 
tions involve simple manipulations of small matrices. 

In Fig. [8] the level scheme for our two models involving 
4 and 6 levels are shown. We first restrict our description 
to the 4-level model as shown in the right part of Fig. [8] In 
this way our description stays compact and matrices are kept 
small. The extension to 6 or more levels follows the same 
scheme. 



1. The 4-level model 

Compared to Fig.J^an additional excited level |3) is added. 
This level corresponds to the rotational level J = and lies 
290MHz below the / = 2 rotational level |1) £J. We work 
in the dressed atom picture and every level \i) is attributed a 
detuning A; (see Fig. [8). Ao is arbitrarily set to 0. The transi- 
tion strengths from the continuum |0) to levels |1) and |3) are 
described by stimulated rates Toi and To3 which are propor- 
tional to the intensity I\ of laser 1 . The transitions between 
the bound levels |2) and 1 1), |3) are characterized by the Rabi 
frequencies £2i2, £223, respectively, which are proportional to 
\fl2- Spontaneous decay from the bound levels leading to 
atomic losses is formally taken into account by introducing 
artificial levels |a,) for each level |/) to which a transition at 
rate ji takes place (not shown in Fig.|8jl. All these couplings 
between different levels are summarized in the symmetric re- 
action matrix K. We arrange the level names in the order 
(0, ai, a2, a3, 1, 2, 3) and use them as row and column in- 
dices. The nonzero matrix elements of the K matrix then read 
K 01 = ^03 = \fT^jl, K iai = K l2 = £2 12 and 

K23 = £^23. Levels 0, ai, a2, a3 are referred to as open chan- 
nels, levels 1, 2 and 3 as closed channels. The reaction matrix 
K is partitioned into open and closed channel blocks, 



K 



K° 



A" oc reads in our case 

/v^W2 v/TV2\ 








K° 







V 



yWJ 



K co is the transposed matrix of K oc , and 



K cc = 



£2i2 
£2i2 £2 2 3 
£2 23 



From K the reduced /T-matrix, 

K led =K 0C (D-K CC ) 



cc\ — l^-co 



is calculated, eliminating the closed channels 1-3, where D 
denotes a diagonal matrix with diagonal elements (Ai , A2, A3). 
This determines the unitary 4x4 scattering matrix S 



red\-l 



S=(l + iK Iea )(l-iK Iea ) 

From the matrix elements Sjj of S the trap loss coefficient Ki ne \ 
is calculated by 



-Kind — 777r5Z I^OaJ 



nti 



l - |Soo| : 



where fi = mRb/2 is the reduced Rb mass and Tiki the relative 
momentum of the colliding atoms. The scattering length is 
obtained from Sqo via 



1 3(S 00 ) 
2ki 91 (S ( 



00 ) 
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where 9t(Soo) an d 3 (Soo) denote the real and imaginary parts 
of Sqq, respectively. 

In the limit of small relative momentum ftki and small cou- 
pling strengths To,- -C Yu Km S \ and the light induced change of 
scattering length a — «b g are independent of k\ because the To, 
are proportional to k\ ( Wigner threshold regime) flfill . 

2. Extension of the 4-level model 

The 4-level model neglects that laser 1 (of which the in- 
tensity is typically five times greater than that of laser 2) also 
couples the levels |2)-|1) and |2)-|3). However, this coupling 
should be taken into account since laser 1 is not far detuned 
from these transitions (see Fig. [HJ due to the small binding en- 
ergy of state |3) (636 MHz x h) which is comparable to typical 
detunings A;. It mainly leads to broadening and light shift- 
ing of level |2). The additional coupling can approximately 



be taken care of by adding another two auxiliary levels |1') 
and |3') with detunings Aj; = A] + A 2 + 2k x 636MHz and 
A 3 / = A3 + A2 + 2k x 636 MHz as shown in Fig. [8] The cou- 
pling strengths Q. ]/ 2 and i2 2 3/ are fixed by = Clny/h/h 

and £2 23 / = ^hzy/h/H- Compared to the 4-level model no 
new fit parameters are introduced. We can calculate Ki ne \ and 
the scattering length a following the same recipe as for the 
4-level model, only with larger matrices. Fitting the data in 
Figs. this extended model produced much better results 
than the 4-level model. For completeness we give here the 
fit parameters which were used in the calculations in Figs.yJ- 
□ (/1 = 300W/cm 2 and I 2 = 60W/cm 2 ): T m /2K = 42kHz, 
r 03 /27r = 8kHz, Q. x2 /2k = 32MHz, £l 23 /27r. = 12MHz, 
Y\/2n = 25 MHz, Yi/2K = 22MHz, y 2 /27T = 2MHz. We used 
k{ = 2.5 x 10 /n -1 . Due to the limitations of ourmodel, these 
fit parameters should not be mistaken as the true values of the 
corresponding physical quantities. 
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